Testable signatures of quantum non-locality in a two-dimensional chiral p-wave superconductor 
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A class of topological excitations - the odd-winding number vortices - in a spinless 2D chiral p-wave (jp x + 
ip y ) superconductor trap Majorana fermion states in the vortex cores. For a dilute gas of such vortices, the 
lowest energy fermionic eigenstates are intrinsically non-local. We predict two testable signatures of this unusual 
quantum non-locality in quasiparticle tunneling experiments. We discuss why the associated teleportation-like 
phenomenon does not imply the violation of causality. 
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Introduction: The intriguing physics of the Majorana 
fermion zero modes at the vortex cores of a spinless p x + ip y - 
wave superconductor/superffuidhas recently attracted consid- 
erable attention dSSHHSSSl- This is because these 
zero modes endow the vortices with non-Abelian braiding 
statistics, and the braiding of such vortices 19J] can, in princi- 
ple, be exploited to build a fault-tolerant topological quantum 
computer (TQC) ifioll. Recently there has been considerable 
evidence II 1 lL 1 1 211 that the symmetry of the superconducting 
order parameter in strontium ruthenate (Sr2Ru04) B13[ |l4ll is 
spin-triplet p x + ip y . It has been proposed in Ref. □ that 
when a thin film of this material is subjected to a magnetic 
field exceeding ~ 200 G, the most energetically favorable vor- 
tices enclose flux hc/Ae rather than hc/2e, and at the core of 
such a half vortex there is a Majorana zero mode. For other 
proposals for the superconducting state in this material, see 
Reft. EES. 

Some heavy fermion materials, such as UPt3, 
are likely to be spin-triplet superconductors as well. More- 
over, with the recent observations of p-wave Feshbach reso- 
nances in spin-polarized 40 K and 6 Li atoms in optical traps 
1 17. Till , these systems are promising candidates for realizing 
p-wave condensates as well [19, 20]. Thus, vortices satisfying 
non-Abelian statistics are tantalizingly close to experimental 
reach. In the following, for brevity, we will use the term "su- 
perconductor" to indicate both superconductor and superfluid. 

An inevitable consequence of the existence of Majorana 
zero modes is that the fermion quasiparticle excitations are 
inherently non-local. The zero energy bound state at the 
vortex core is described by the Majorana fermion operator 
7^ = 7 [!]■ Other bound states are separated by an energy 

A 2 

gap u>o ~ -j 11 - <C An, where An gives the gap in the bulk 
and ep is the Fermi energy. If the system has two Majorana 
zero modes, 71 and 72, one can define a composite operator 
= (71 + 172) /V2, and its hermitian conjugate q, which 
satisfy Fermi anticommutation relations. If an electron is in- 
jected into the system with an energy <C (Jq, it can only go 
to the excited state described by q and q*. This intrinsic non- 
locality of the quasiparticle wavefunction was first discussed 
in Refs. l2il.l22l for a ID quantum wire embedded in a p-wave 
superconductor. In this paper, we use it to make interesting 
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FIG. 1: (a) Two quantum dots in tunneling contact with two vortices, 
(b) A quantum dot and a STM tip in tunneling contact with two vor- 
tices. 



predictions for quasiparticle tunneling experiments in a 2D 
Px + ipy superconductor. 

To distinguish the non-locality discussed above from that 
in a quantum double well problem, let \l) and |r) represent 
the two states localized in the left and the right well. In 
the presence of quantum tunneling, ~t(cjc r + h.c), the 
ground and the excited states are given by the symmetric 



and the anti-symmetric states, 
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r) |0> and 
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= (cj — cj) jO), respectively. Thus, it might appear that 
in this case also the quasiparticle operators are delocalized 
between the two wells even when they are very far separated. 
However, if we inject an electron into the left well, it will 
enter into the linear combination of Is) and la), the state lo- 



calized in the left well 112 U 12211 . Only after a time of the order 
of h divided by the tunnel splitting does the electron acquire 
appreciable probability of appearing at the right well. Further, 
because of the exponentially small tunnel splitting, by adding 
a small perturbation, e.g., clc l cl, to the Hamiltonian, the 
eigen wavefunction will become localized near the left or the 
right well. Motivated by these considerations we consider the 
following thought experiments to test the explicit non-locality 
of a quasiparticle excitation in the Majorana fermion system. 

Teleportation: We consider the set up of Fig. 1(a) where 
two voltage biased quantum dots are in tunneling contact with 



2 



two non-Abelian vortices of a p x + ip y superconductor. The 
effective Hamiltonian describing the system at energy < < ojq 
is given by, 

H = ^ e a1pil/ja + iJjLjR ~ V%t L ~{ L (%p\ - 1p L ) 
a=L,R 

-V2t RlR {^ R -^ R ) (1) 

Here, jl and are the two Majorana fermion operators 
and ipL and i/jr are the electron annihilation operators for 
the lowest energy state of the quantum dots. The parameter 
J describes the hybridization between the Majorana fermion 
modes. It is exponentially small if the distance between the 
vortices is large. 

After introducing fermion operators = (7^ + ijl)/ V%, 
q = (jR ~ «7l)/v / 2, and performing the local gauge trans- 
formation ipL — > itpL, Eq. (Q]) becomes [23] 

H = ^ e Q i/)^ Q + Jq^q + ^ t a (q^ ifj a + h.c.) 

a=L,R a=L,R 

-t L {*>py + h.c.) + i fl (V>V + h.c). (2) 

In this form, it is clear that the tunneling terms only conserve 
charge modulo 2e. (The extra or deficit charge is taken away 
or provided by the superconducting condensate.) The low- 
est energy Hilbert space upon which Eq. (f2j) acts is eight di- 
mensional. Its basis states are labeled by n q , Ur), where 
n L,q,R = 0,1 denote the number of particles in the single par- 
ticle states associated with the t/>l, g and ip R operators, respec- 
tively. Because of charge conservation modulo 2e, this eight 
dimensional Hilbert space decouples into two four dimen- 
sional ones spanned by {|1, 0,0), |0, 1,0), |0, 0, 1), |1, 1, 1)} 
and {|0, 1, 1), 1 1,0, 1), 1 1, 1, 0), 1 0,0, 0)}, respectively. In the 
demonstration of the teleportation phenomenon, we shall con- 
sider one extra electron at the left dot at time zero and see how 
it is transported to the right dot at a later time. For that pur- 
pose, we only need to consider the subspace where the number 
of electrons is odd. In this subspace the Hamiltonian is a 4 x 4 
matrix given by 



H = 
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To simplify the calculations, in the following we shall con- 
sider the situation where J — > + , = t R = t, e R = 0, and 
study the eigenstates of Eq. (HJ as a function of . Physically, 
we expect that when e_L is large and negative, the extra elec- 
tron initially at the left dot will stay there. On the contrary, by 
tuning £l to zero, the extra electron can delocalize by tunnel- 
ing to the right dot via the conduit provided by the extended 
q-state. Notice that this is very different from ordinary quan- 
tum tunneling of the electron which is non-zero only when J 
is appreciable. 

A simple calculation shows that if the initial state is 
1 1, 0, 0), the probability of observing an electron at the left 



and the right dots at a later time T is given by 

- ,;ffT |l,0,0)| 2 + |(l,l,l|e-^ T |l,0,0)| 2 
9t 2 

= 1 - 



P L {T) = |(l,0,0|e- 
2t 2 



At 2 



1 - cos (TJM 2 + e|) 



P R {T) = |<0,0, l|e-^ J |l,0,0)| 2 + |(1,1,1| 
= sm 2 (T<). 



-iHT 1 



1,0,0)| 2 
(4) 



Eq. dUl predicts that Pl — * 1 for large negative e^. Therefore, 
as expected, the electron initially on the left dot will remain 
there. What is more important is the fact that the probability 
of finding an electron on the right quantum dot is independent 
ofe L \ 

If one separates P R {T) into the sum of the following prob- 
abilities, 



P m (T) = |(0,0,l| e -^ T |l,0,0)| 2 

4ffT |i,o,o)| 2 , 



P R2 {T) = |(1, l,l|e 
one would obtain, 



(5) 



Prx{T) 



it 2 



3in 2 (7Y) sin 2 ( 



T^/At 2 + e 



Pri(T) = sm 2 {Tt)-P R1 (T). 



(6) 



The first term, P R \(T), is the probability of transporting the 
electron from the left to the right dot without disturbing the su- 
perconductor. For €l = it is equal to P R i{T) — sin 4 (Tt), 
hence it equals to unity when T = it /2t. The fact that the 
above time is finite even in the limit of J — * + suggests that 
it is independent of the (large) separation between the vortices, 
thus justifying the word "teleportation". 

Note that the total probability of observing an electron at 
the right dot, Eq. (0), is independent of what's being done at 
the left dot. Since there is no way to distinguish between the 
processes responsible for P R i and P R 2 by a measurement 
only on the right dot, there is no way to know if the electron 
has been teleported from the left dot or has arisen out of 
the condensate. As a result, no information can be sent 
from the left to the right dot, hence causality is maintained. 
To prove the existence of teleportation, it is necessary to 
differentiate the processes responsible for P R i and P R 2- For 
example, one could perform a coincident measurement of 
both the occupation numbers of the left and the right dots 
at time t. This measurement will give zero for the process 
responsible for P R \ and non-zero for that responsible for 
P R 2- However, since a 'classical' exchange of information 
(the result of the coincident measurement) is necessary, there 
is no superluminal transfer of information in the observation 
of the teleportation effect. This experiment is possible in 
principle, but will be difficult in practice. In the following, 
we describe a much easier STM experiment which reveals the 
non-local nature of the q-quasiparticle state discussed earlier. 

Action over distance: Now let us consider the set up in 
Fig. 1(b) where a voltage-biased quantum dot is in tunneling 
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contact with the left vortex, while an STM tip probes the dif- 
ferential conductance dl/dV of the right dot. The purpose is to 
study the effect of biasing the left dot on the tunneling curve 
at the right dot. Again one might expect a non-zero effect due 
to the extended q-state. 

In computing the STM spectral function, we use the formu- 
lae, 

A(u) = J2\(a-n\^ R K)\ 2 S(uj - + E%) (7) 

a 

for to > 0, and 

A{u) = J2\(^HrK)\ 2 S(u; - + E%) (8) 

a 

for oj < 0. Here rj = ±1 labels the even and odd particle 
number subspaces respectively, and \a n ),EH are the eigen- 
states and the eigenenergies of Eq. © (a = indicates the 
ground state). In the explicit calculation, we diagonalize the 
following 8x8 Hamiltonian matrix, 



H 



E x 
H 2 



(9) 



where E\ is the 4x4 matrix given in Eq. (f3]) and H 2 is given 

by, 
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In this 8— dimensional space ip a is represented by 
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For J = 0, straightforward computation gives 

i r[2t 2 + r 2 + M(M-2t)] 



A{w) 



2ir (lu 2 + T 2 )[{\uj\~2t) 2 +T 2 } 



(13) 



where T is a Lorentzian broadening of the energy levels. The 
fact that the result is independent of el is expected because, 
for zero J, there is no communication between the two Ma- 
jorana fermion states associated with an d Tfl- As a result, 
Eq. (Q]l decouples into two independent quantum dots in tun- 
neling contact with two independent vortices. Because of this, 
the tunneling spectrum of the right dot is independent of the 
left dot. 

Next let us consider the case of non-zero J. In this case, 
the expression for A(u>) is lengthy and we shall just present 
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FIG. 2: The tunneling curves for J = 0.5t. Dashed curve: = 0, 
solid curve: e L =t. 



the numerical results for J — 0.5t. In Fig. 2 the dashed curve 
is the tunneling spectrum for el = while the solid curve 
is that for el = t. In constructing this figure we used V — 
0.15t. It is clear that the tunneling spectrum at the right dot 
can be qualitatively altered by tuning cl only on the left dot. 
The most surprising result is that, when el ^ 0, there is a 
discontinuity in A(u>) at u> = 0! The presence of zero energy 
spectral weight is due to the relation: 



E~ v 



E n 



The discontinuity arises from the fact that 

|<o_,|^|o,>|VKO-^b|o,>| s 



(14) 



(15) 



Eq. ( TPfl ) is due to the fact that, for er — 0, H\ and H 2 have 
identical eigenvalues. As for Eq. ( fT5l ). it can be shown analyt- 
ically that 

KO-.I^Jo^p-Ko^l^jo,)! 2 

= |(0 I! ^fl|0_,)| 2 -|(0_ I! ^fl|0,)| a 
Je l 



a/16* 4 + e 2 L (At 2 + J 2 ) 



(16) 



which yields zero (no discontinuity) if either J or ej, is zero. 

A discontinuity in the tunneling spectrum is very unusual. 
The presence of a non-zero tunneling spectral weight at to = 
requires the absence of an energy gap. This occurs, for exam- 
ple, in metals and Anderson insulators. The discontinuity in 
the spectral weight at uj = then requires |(0jv-i|V'R|0iv)| 
|(0Ar|^jj|0jv+i)|, where \0m) is me M-particle ground state. 
This would not occur in a metal or an Anderson insulator. 

Experiments in atomic superfluid: The schemes proposed 
above for testing the non-local properties of the Majorana 
zero modes can be implemented in a (p x + ip y )-wave atomic 
superfluid through a recently proposed laser probing scheme 
B24l 12511 . This scheme is similar to the tunneling microscopy 
of a superconductor in that it relies on induced tunneling be- 
tween a superfluid and a normal phase. Here, we assume that 
the atoms in the superfluid are prepared in a hyperfine state 
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1 1), while the atoms in the other hyperfine state | f ) do not par- 
ticipate in the p-wave pairing, and therefore constitute the nor- 
mal phase. The interference between the superfiuid and nor- 
mal phases can be realized through a two-photon Raman pro- 
cess that couples the two states |J.) and |f) of atoms with two 
local laser fields. An optical dipole trap with spin-dependent 
potential depth can be used to hold the normal atoms near the 
vortex cores. The bias voltage el can be adjusted by vary- 
ing the intensity of the optical dipole trap, which changes the 
energy of the normal atoms with respect to atoms in the super- 
fluid. The number of atoms in the normal phase can be mea- 
sured through fluorescence signals of resonance lasers. The 
tunneling strengths ti, and t R can be adjusted by varying the 
intensities of the Raman lasers fields, which changes the cou- 
pling strength between the two hyperfine states. In the exper- 
iment, to detect the action over distance, the spectral function 
A {lo) can be measured by varying the detuning of the laser 
fields to the atomic transition 5 = ujl—lqa 11241 125II . where lo l 
is the frequency difference between the two laser fields for the 
two-photon Raman process, and loa is the frequency splitting 
between the two hyperfine levels \\) and ||). In this way, the 
proposed experiments to test the quantum teleportation and 
the action over distance can be performed on an atomic su- 
perfiuid as well. The mathematical illustrations, including the 
form of the Hamiltonian and the measurable quantities, re- 
main the same as in an electronic superconductor. 

Experimental feasibility : Here we comment on the cut-off 
energy scales below which the interesting effects discussed 
above would be observable. The tunneling amplitudes and 
t R (the bias voltages) need to be much smaller than Uq. For 
Sr2Ru04 (assuming it is a spin-triplet p-wave superconduc- 
tor) using T c ~ 1.5 K and the effective fermi temperature 
T F ~ 50 K iQjl], this requires t R , t L < u ~ 50 mK. In the 
case of Feshbach superfluids, where Ao ~ cf in the unitary 
regime [26], ujq ~ ep and the effects due to the Majorana 
fermions should be observable below the scale of the Fermi 
temperature itself. It is also important that the experiment be 
done on a time scale short compared to the decoherence time. 
The relevant processes are (1) thermal excitation of a fermion 
zero mode to a higher-energy state, and (2) quantum tunnel- 
ing of the zero mode from the vortex probed in our experi- 
ment to some other vortex induced by disorder. Process (1) 



merely need T <C £f- The second process can be avoided 
by performing the experiment with vortices which are well- 
isolated from any other vortices (which can be detected by 
various imaging techniques). 

To conclude, the observation of our predicted nonlocal 
quantum entanglement behavior would be direct evidence 
supporting the existence of non-Abelian topological anyons 
in 2D chiral p-wave superfluids and superconductors. 
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will occur with probability 



-uio/T 



For T < 50mK, this 



will be small in Sr2Ru04 (note that disorder may also sup- 
press the excitation gap ujq locally, which would require the 
temperature to be reduced further); in a cold atom system, we 
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